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ABSTRACT 

We present compelling evidence for the complexity of the Fornax dwarf spheroidal. By disentan- 
gling three different stellar sub-populations in its red giant branch, we are able to study in detail 
the dependence between kinematics and metallicity. A well-defined ordering in velocity dispersion, 
spatial concentration, and metallicity is evident in the subpopulations. We also present evidence for 
a significant misalignment between the angular momentum vectors of the old and intermediate-age 
populations. According to the HST measurement of Fornax's proper motion, this corresponds to 
counter-rotation. These ingredients are used to construct a novel evolutionary history of the Fornax 
dwarf spheroidal, characterized as a late merger of a bound pair. 

Subject headings: galaxies: kinematics and dynamics — Local Group — dwarf — individual: Fornax 
dwarf Spheroidal 



1. INTRODUCTION 

As the size and quality of datasets on the local pop- 
ulation of dwarf Spheroidals (dSphs) increase, evidence 
for the complexity of these systems is intensifying. The 
observed variety of properties and structures gives each 
dwarf its own distinctive character. Their star formation 
histories, for instance, are well known to differ signifi- 
cantly. These often give rise to multiple generations of 
stars, but the mechanisms that drive an intermittent star 
formation must be investigated on a case-by-case basis. 
Being able to untangle the properties of such stellar sub- 
populations is crucial in advancing our understanding of 
how they came in place. 

After the Sagittarius dwarf, Fornax is the second most 
luminous dSph orbiting the Milky Way and the only 
one (in addition to Sagitt arius) to posses an as sociated 
Globular Cluster system ()van den Berghl Il998ft . This 
is not Fornax's only peculiarity: for example, signifi- 
ca nt asymmetrie s in the isophotes were recorded as early 
as jHodgd (11961D. lEskridgd (119881) and then confirmed 
llrwi 



in llrwin fc Hat zidimitriou ( 1995 Q . Recent deep photo- 
metric surveys dStetson et al.l 11998; Sav iane et al.l 120001 
iBattaglia et al.ll2006l B06 in the following) have revealed 
a rich and prolonged star formation, comprising old stars 
(> 10 Gyr), a dominant population of intermediate age 
stars, as well as stars of only a few hundred Myr. 

All s ystematic studies of Fornax's star formation his- 
tory (IColeman fc de Jonel [20081: iDel Pino et ail [20lTI: 
Ide Boer et al. I 12012ft record a significant starburst at 
approximately 4 Gyr ago. The trigger of this activ- 
ity, though, is still debated. On the one hand, at least 
three stellar overdensities have been identified in F ornax 
(|Coleman et al.l l200l [20051 Ide Boer et aT~l l2012t ). and 
interpreted as remini scent of the shell featu res observed 
in elliptical galaxies (|Malin fc Carterlll980l) . In this pic- 
ture, Fornax collided with a low-mass sub-halo in the 
relatively recent past and swallowed it. On the other 
hand, the details are not free from difficulties, especially 
caused by the inferred metallicit y and age of the stars 
composing the innermost clump {Olsz ewski et all 120061 : 
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IColeman fc de Jongil2008ft . These are about 1.5 Gyr old, 
which makes the timing of the global starburst problem- 
atic. They are also relatively metal-rich, which in turn 
does not implicate an external origin for their gas. Fur- 
thermore, a recent collision is quite unlikely when the 
energetics of encounters in a yiriali zed Milky Way halo 

is considered (|De Riicke et al.ll200l). 

Afte r the pioneering work of IWalker fc Pefiarrubial 
(| 2 1 lh (WPll in the following), in this Letter, we pro- 
vide new evidence for Fornax's complexity by disentan- 
gling different stellar subpopulations in the red giant 
branch (RGB) and by studying in detail their kinemat- 
ics. Rather than a division into two subpopulations 
()Saviane et al.ll2000L B06, WPll), one into three is pre- 
ferred by the data. The identified subpopulations show 
significant differences in line-of-sight (LOS) velocity dis- 
persion and higher order moments. Even more striking 
is the incompatibility between their detected rotations, 
with the relative rotation vector inclined ~ 40° with re- 
spect to the isophotal major axis. Finally, we use these 
new ingredients to inform a novel formation scenario for 
the Fornax dSph. 

2. THREE BAYESIAN SUBPOPULATIONS 

A step forward in modelling the superposition of 
chemodynamically distinct stellar populations in dSphs 
has been achieved very recently by WPll. They de- 
vise a Bayesian technique that is able to quantify and 
separate the distributions of projected radii, velocities 
and metallicities of distinct subcomponents. Here, we 
use this technique extensively. Given a spectroscopic 
dataset, each star is associated to a set of probabilities 
of membership, one to each of the identified stellar sub- 
populations. This is useful, as it allows us to avoid any 
rigid color or metallicity cut, which inevitably degrades 
the data. Instead, these probabilities can be used in de- 
riving robus t and detailed kinematics for each subcom- 
ponent (e.g.. lAmorisco fc Evans 2012, A12 in the follow- 
ing). Each population is assumed to have a Plummer 
surface density, a Gaussian metallicity distribution and 
a Gaussian LOS velocity distribution (the latter hypoth- 
esis is relaxed in Section 4). Also, this technique is able 



2 



to correct for the selection effects introduced by the spa- 
tial distribution of the spectroscopic sample, hence all 
inferred quantities - such as half-light radii and fractions 
of giants - are genuine global quantities. 

WP11 apply this technique to the ( approximately 
2500) red giants in the lWalker et al.l (2009) spectroscopic 
dataset (W09 in the following), and look for a decom- 
position in, at most, two subcomponents. Indeed, they 
distinguish a dominant, colder and more concentrated 
metal-rich (MR) stellar population, from a more diffused, 
kinematically hotter metal-poor (MP) population. How- 
ever, we reckon that such a division into two subpopu- 
lations is not stable with respect to removal of the kine- 
matic information. If only the spatial distribution and 
the metallicities of the giants are used (i.e. discarding the 
LOS velocities), a different division is recovered. This is 
not the case if the presence of a third population is al- 
lowed, as shown in Table 1. In fact, a three-population 
division is a significantly better description: since the 
two models are nested, the comparison is straightforward 
and a simple m aximum- likel i hoo d ratio (MLR) is suffi- 
cient. Following lEadie et al.l (|1971l ). we find that despite 
the penalty given by the additional free parameters, the 
probability of obtaining such a MLR by pure chance is 
less than 3 x 10~ 6 . Additionally, no evidence is found in 
the data to justify the addition of a fourth population. 

While the clear correlations between metallicity, spa- 
tial distribution and velocity dispersion are naturally 
maintained in passing from two to three subpopulations, 
the data prefer to accommodate separately a small frac- 
tion of MR giants, which in turn have colder kinemat- 
ics and are more centrally concentrated. A division into 
three subpopulations also best traces Fornax's star for- 
mation history, allowing at the same time for stars that 
are older and younger than the bulk of intcrmcdiatc- 
metallicity (IM) giants. In this sense, the presence of 
such a feature can be easily understood, as the abundant 
1-2 Gyr old population of stars, traced for example by 
the blue loop stars (B06 or de Boer et al. 2012), are also 
expected to sneak into the RGB. 

The top panel of Fig. 1 shows the final metallicity sub- 
division of the W09 dataset. The MR peak is scarcely vis- 
ible against the strong and broad IM population, whilst 
the middle panel shows the cumulative mass distribu- 
tions of the three populations, illustrating the impor- 
tance of combining evidence from the metallicity and 
the spatial distribution of the giants. A confirmation 
of the statistical importance of the MR feature comes 
from the independent spectroscopic dataset of B06. This 
has a different spatial distribution and contains at least 
400 giants that do not appear in the W09 dataset. Fur- 
thermore, metallicities have been calibrated from spectra 
in the Ca II triplet region, and are affected by smaller 
uncertainties, resulting in a significantly less broadened 
metallicity distribution in the lower panel of Fig. 1. We 
apply the same statistical technique and again obtain 
that a division into three populations is a significantly 
better description (the probability of obtaining such a 
MLR by chance is < 5 x 10~ 4 ). The MR feature sits at 
[Fe/H] fx —0.65, and is more concentrated than the other 
two subpopulationtQ. 



1 We do not record information on the fractions or scale radii of 




0.3 0.4 0.5 

Z Mg /A 




R/kpc 




Fig. 1. — Upper and middle panels: metallicity distributions 
and cumulat i ve ma ss distributions (one sigma errors) for the 
IWalker et al. (2009) dataset . Low er panel: metallicity distribu- 
tions for the Battaelia et al. (2006) dataset. Color-coding is asso- 
ciated with increasing metallicity. 
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Fig. 2. — Top-le ft: rotation pattern from the datasets of 
IWalker et alj l|2009l ). one and two sigma contours for each stel- 
lar pop ulation, color-coding as in Fig. 1; in black, Fornax's PM 
as from Piatck et al. (2007). Top-right: rotation for inner (dot- 
ted ellipses) and outer (dashed ellipses) stars for the IM and MP 
population. Lower panel: the intrinsic rotation signal (after sub- 
traction of the PM) along the di rection of the re l ative rotation, 
as a function of metallicity, for the Battaglia et al. (2006) dataset. 
Datapoints of growing sizes (and darker shades) are obtained from 
larger subsamples. 



the three subpopulations as obtained from the B06 dataset. Unlike 
the metallicity, these are not reliable unless the spatial sampling 
bias is taken into account. 
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TABLE 1 

TWO AND THREE POPULATION DIVISIONS OF THE FORNAX DWARF 



Stellar 
Population 


[in A\ 


StD(SMg) 
[in A] 


Rh 
[in pc] 


Fraction of 
Population 


(«) 
[in kms 


— "1 1 

• ] 


Metal-poor 


0.31 ± 0.025 


0.09 ±0.015 


860 


± 


60 


0.40 


± 


0.09 


14.4 


± 


0.6 




0.24 ±0.015 


0.04 ±0.015 


105C 




:90 


0.19 


± 


0.04 








Metal-rich 


0.47 ±0.01 


0.081 ± 0.005 


556 


± 


24 


0.60 


± 


0.09 


10.0 


± 


0.6 




0.447 ±0.007 


0.091 ± 0.004 


598 


± 


16 


0.81 


± 


0.04 








Metal-poor 


0.27 ±0.018 


0.07 ±0.01 


888 


± 


50 


0.31 


± 


0.06 


14.3 


± 


0.6 




0.26 ±0.015 


0.05 ±0.01 


935 


± 


65 


0.27 


± 


0.09 








Intermediate 


0.46 ± 0.01 


0.058 ± 0.006 


605 


± 


30 


0.56 


± 


0.05 


11.3 


± 


0.7 




0.453 ± 0.009 


0.065 ± 0.008 


610 


± 


25 


0.63 


± 


0.07 








Metal-rich 


0.53 ±0.035 


0.11 ±0.02 


480 


± 


45 


0.13 


± 


0.04 


8.6 


± 


1 




0.56 ± 0.05 


0.10 ±0.025 


437 


± 


55 


0.10 


± 


0.04 









Note. — Top-part: detailed results for a two-population division. Lower part: a third population is allowed. Results are collected for 
the case in which the kinematic information is retained (first line for each population), or discarded (second line). 



3. A COMPLEX ROTATION PATTERN 

Since Fornax is an extended object, those stars that 
are misaligned with the barycenter have their LOS veloc- 
ities affected by a geo metrical projectio n of the systemic 
proper motion (PM) (|Feast et alj|196l| ). This projection 
is equivalent to an apparent solid-body rotation; hence, 
any rotation signal detected in the LOS velocities is ac- 
tually the superposition of any intrinsic and apparent 
rotation: 

^tot = ^int + ^app , (1) 

where f2 app is proportional to the systemic PM /2 (e.g. 
eqn (3) in A12). For clarity, our frequencies are defined 
so that 

Vlos(£Maj, £Min) = ( ) " (^Maj , ^Min) , (2) 

with all cartesian components referring to the major and 
minor axes of the isop hotes . 

I Walker et al.l (|2008| ) have recently measured the total 
rotation signal in Fornax. They find it in agreement with 
the direct as trometric det e rmina tion of Fornax's PM, as 
obtained by iPiatek et all (|2007l henceforth P07) using 
HST data. It was then natural to interpret this result as 
evidence that Fornax has very little intrinsic rotation in 
the LOS direction. Indeed, both the W09 and the B06 
datasets are consistent with this interpretation, although 
this simple picture remains valid only as long as all the 
giants in each dataset are considered together. When 
we draw subsamples, we detect a strong dependence of 
the total rotation signal on metallicity, with a richness of 
detail that is surprising for such a small system. 

In order to measure total LOS rotations Ototj we 
adopt the same m aximum-likelihood technique as in 
iWalker et al.1 (|2008| ). Our probabilities of membership 
allow us to measure the rotational properties of each 
subcomponent. As a compromise between size of the 
subsamples and excessive contamination, we only con- 
sider giants with at least 50% chance of belonging to 
a specified population, but our results are qualitatively 
stable against changes in this choice. Results are shown 
in Fig. 2. 

The measure we obtain for the MR feature has a large 
uncertainty, due to the limited number of stars (~ 80 
giants with pmr > 0.5), as well their spatial concentra- 
tion. On the other hand, the rotations of the MP and 
the IM populations are precise enough to clearly detect 



a relative rotation: 

fUi - ^ - = i$J - . (3) 

This is inclined at ~ 40° with respect to the major 
axis and provides a natural justification to the angu- 
lar behaviour of the asymmetries detected in the LOS 
velocity distribution by A12 (see their Fig. 12). Its 
magnitude is independent of the systemic PM and is 
non- negligible, amounting to at least 3 kms _1 kpc _1 - 
as a comparison, Fornax's ti dal radius is ~ 2.85 kpc 
(jlrwin fc H atzidimitriovj fl995t l. Furthermore, if we as- 
sume that the PM measured by P07 is correct, then this 
relative rotation is also a counter-rotation. 

The B06 dataset confirms this picture. In fact, since 
their metallicites are more accurate, it provides an even 
more detailed description, as shown in the lower panel 
of Fig. 2. A comparison with Fig. 1 shows that peaks 
at positive and negative values of fl^o trace the peaks 
in the metallicity distribution, respectively, of the MP 
and IM populations. Nonetheless, additional structure is 
also visible. The most significant lies in the MP end of 
the distribution, where the rotational velocity SX/^° drops 
again at negative values. This feature is in fact rather 
thought-provoking because it cannot be interpreted as 
an effect of contamination of the MP population from 
the other subpopulations. Still, it also appears in the 
W09 dataset. In fact, additional indications that the 
MP population may indeed have a richer structure come 
from the top-right panel of Fig. 2, in which the W09 
subsamples pertaining to the IM and MP populations 
are split into two circular annuli each. While for the IM 
the innermost giants (dashed ellipses) share the same ro- 
tation as the outermost ones (dotted ellipses), we detect 
a significantly more pronounced relative rotation in the 
innermost MP giants (approximately within 1 kpc). 

4. DISPERSIONS AND HIGHER ORDER MOMENTS 

A12 have recently developed a Bayesian implementa- 
tion for the extraction of the shape information of line 
profiles from discrete data. We apply that here and, for 
the first time, measure resolved kinematics of the stellar 
subpopulations in Fornax. It is probably worth mention- 
ing that, in order to obtain reliable information on the 
higher order moments, it is necessary to use the results 
on the division into three populations as obtained by dis- 
carding the kinematic information (see Section 2). 

The left panel of Figure 3 collects the velocity disper- 
sions profiles - in circular annuli - for the three iden- 
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Fig. 3. — Velocity dispersion and symmetric deviations in circular 
annuli for the three stellar subpopulations. 

tilled subpopulations, with the same color-coding as in 
Fig. 1 (additional red-dotted points use a smaller bin- 
size for the MR population). It is reassuring to find 
that a Bayesian dissection only based on spatial posi- 
tions and metallicities also yields significantly different 
kinematic profiles. In fact, also the shape of the ve- 
locity distributions shows significant differences. We re- 
call that the symmetry parameter s is analogous to the 
Gauss- Hermite coefficient /14, and quantifies - with re- 
spect to a Gaussian (s = 0) - how 'peaked' (s > 0) 
or 'flat-topped' (s < 0) the velocity distributions is (see 
A12 for further details). A12 recorded a preference for 
a mild flat-toppedness in Fornax (s w —0.3), which we 
can now interpret as being due, at least in part, to the 
detected relative rotation. Instead, once isolated, the IM 
population prefers quasi-Gaussian profiles, with a mild 
tendency for positive values of s. 

5. A MERGER OF A BOUND PAIR? 

We have presented a detailed view of the rich kine- 
matical properties of the apparently puny Fornax dSph. 
Both spectroscopic datasets (B06 and W09) support a 
Bayesian division into three populations of its RGB, with 
a new MR feature being identified at [Fe/H] m —0.65, 
containing approximately a tenth of the giants. Both 
datasets also show a strong dependence between stellar 
kinematics and metallicity. This includes a well-defined 
ordering of the three populations in velocity dispersion 
and spatial concentration, but it is most dramatic in a 
clear misalignment between the angular momentum vec- 
tors of the MP and the IM populations. The magnitude 
of this difference is independent of the systematic PM of 
Fornax, but, if we assume that the measurement by P07 
is correct, then this misalignment also translates into a 
counter-rotation. The populations are then rotating in 
opposite directions along an axis inclined at ~ 40° with 
respect to the major axis. 

Together with other evidence in the literature, the 
abundance of features complicates the task of a co- 
herent reconstruction of Fornax's evolutionary history. 
For instance, a merging scenario seems to be natu- 
rally inv oked by the shell features p resent in the pho- 
tometry (jColeman et all [20041 . [2005). Even so, this is 
made less convincing by the realization that the gas 
fueling the stars in such overdensities is more likely 
from Fornax itself, given its significant pre-enric hment 
([Olszewski et all 120061 : fColeman fe de Jongjl2008D . Fur- 
thermore, the merger scenario has to overcome the diffi- 
culty of having a collision between two sub-haloes orbit- 



ing t he Milky Way at such late times (|De Riicke et alj 
l2004f ). On the other hand, it is almost impossible to rec- 
oncile the misalignment between the angular momenta 
of the MP and IM population without interactions. If 
we are skeptical about the PM measurement by P07, we 
might relax the argument for counter-rotation, and then 
invoke a tidal origin for the different rotational properties 
of the MP stellar population. Nonetheless, the fact that 
it is actually the innermost MP stars that have a faster 
relative rotation rather than the outermost ones argues 
against this, since those stars sit at the same distance 
from the center as most of the IM stars. 

Furthermore, supporting evidence is provided by the 
fact that the spatial distribution of the youngest stars 
in Fornax seems to have knowledge of th e direction of 
the relative rotation. IStetson et al.l (|1998D note that the 
main sequence stars "display a clearly flattened distri- 
bution on the sky, with a long axis in the east-west di- 
rection" (see their Fig. 12), hence in perfect alignment 

with O r ei. It is then appealing to advocate that this cen- 
tral structure, composed by stars that are younger than 
a few 100 Myr and whose geometry is uncorrelated with 
Fornax's isophotes, is reminiscent of an elongated bar. If 
so, its spatial alignment is also the direction of its intrin- 
sic rotational velocity. In turn, this is aligned with the 
intrinsic rotations of the MP and IM populations only if 
the PM is not significantly different from P07. 

We suggest that a solution for these inconsistences is 
that Fornax is a merger of a bound pair. First, if Fornax 
and its perturber have some relative orbital motion at 
early times, say with frequency f2 rb, then interactions 
with the old MP population provides an efficient means 
of heating a nd preferentially stripping the progradel y ro- 
tating stars p'Onghia et al.ll20ldlRead et al.ll2006l) . As 
a result, whatever the initial internal orbital structure, 
the MP population is left with a counter-rotating signa- 
ture (with respect to f2 rt>)- Second, the nearby pres- 
ence of a companion increases significantly the chances 
of retaining the gas expelled by the first generations of 
stars in Fornax. Its escape is made more difficult by the 
increased total mass and by the gas-dynamics of the in- 
teractions with the companion. As dynamical friction 
eventually causes the decay of the perturber 's orbit, the 
infalling gas naturally has both a high metallicity and 
a prograde rotational signature, hence imprinted in the 
younger stellar populations. This proce ss may be pro- 
longe d, and the localized overdensities (jColeman et al.1 
[2004D and the central bar of main sequence stars may 
represent the last stages before the supply of gas is ex- 
tinguished. 

How massive might such a putative companion be? By 
assuming that the luminosity in the RGB is about one 
half of the total in Fornax, that the luminosity func- 
tions of the three sub-populations are the same, and that 
the metal poorest tail contains about one third of the 
MP giants, we get L m 7 x 10 5 Lq, roughly compara- 
ble to Canes Venatici I, Draco, Sextans, Ursa Minor, 
although this may just be a lower limit. These dSphs 
have masses of ~ 4 x 10 7 within 1.7 half light radii 
(e.g.. lAmorisco fe E vans! [20 111 ), therefore probably more 
than a few 10 8 M e in total. While this is still too low 
to make a late collision probable, it is instead enough to 
significantly perturb the progradely rotating stars in For- 



5 



nax, and to ensure an infall timescale of a fraction of the 
Hubble time. We conclude that this picture collects the 
benefits of the merging scenario while overcoming most 
of its limits. 
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